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1.  Introduction 


The  ultra-wideband  (UWB)  radar  technology  has  many  applications  to  sensing 
scenarios  where  high  resolution  and  typically  low  microwave  frequencies  are 
required,  such  as  sensing  through  the  wall,  ground,  or  foliage  penetrating  radar.  The 
US  Army  Research  Laboratory  has  investigated  this  technology  for  more  than  2 
decades,  with  an  emphasis  on  radar  imaging  applications. 

A  previous  study  by  this  author1  suggested  that  UWB  waveforms  are  not  the  best 
choice  for  Doppler  radar  processing  due  to  the  issue  of  target  range  migration 
during  the  integration  interval.  However,  more  recent  investigations2,3  have 
demonstrated  that  this  problem  can  be  overcome  by  the  correct  choice  of  the 
matched  filter  in  the  case  of  UWB  waveforms.  While  the  ideas,  mathematical 
derivations,  and  numerical  examples  in  this  note  were  developed  independently, 
they  are  very  similar  to  the  results  already  published  by  Le  Chevalier2  and  He  et 
al.3  Although  not  entirely  original,  this  note  was  written  mainly  as  an  update  to  the 
conclusions  of  our  previous  work.1  Additionally,  we  include  a  discussion  of  the 
computational  aspects  of  this  technique  as  well  as  a  detailed  account  of  the 
wideband  ambiguity  function  (AF)  and  related  Doppler  ambiguity  issues,  which,  to 
our  knowledge,  have  not  been  published  elsewhere. 

This  report  is  organized  as  follows.  In  Section  2,  we  develop  the  theory  behind  the 
new  Doppler  processing  technique  using  UWB  waveforms.  Section  3  presents  a 
simple  numerical  example  and  compares  the  new  method  with  the  conventional 
Doppler  processing  approach.  In  Section  4,  we  discuss  both  the  computational 
aspects  of  implementing  the  wideband  matched  filter  and  the  issue  of  Doppler 
ambiguities  in  the  context  of  the  new  processing  scheme.  We  finalize  with 
conclusions  in  Section  5. 

2.  Theoretical  Description 

The  following  symbols  are  used  in  this  and  the  following  sections: 

p  -  complex  envelope  of  the  transmitted  waveform  in  the  time  domain 
P  -  transmitted  waveform  envelope  spectrum  (frequency  domain) 
x  -  transmitted  waveform  (complex  form) 
y  -  received  waveform  (complex  form) 
z  -  demodulated  received  waveform  (complex  form) 
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s  -  2-D  received  signal  in  the  fast  time  domain 
S  -  2-D  received  signal  in  the  frequency  domain 
Q  -  phase  factor  in  the  2-D  point  target  response 
RDM  -  2-D  range-Doppler  map 
t  -  time  (in  general) 
tf-  fast  time 
ts  -  slow  time 

Tr  -  pulse  repetition  interval 
Tcpi  -  coherent  processing  interval 
r-  received  waveform  delay 
r  -  target  range 

R  -  target  range  in  the  middle  of  the  coherent  processing  interval  (CPI) 

/-  RF  frequency  (Fourier  pair  of  fast  time) 

/o  -  carrier  frequency 
fd  -  Doppler  frequency 
B  -  pulse  bandwidth 
v  -  radial  velocity  of  the  target 
c  -  speed  of  light 

p  -  variable  representing  the  Fourier  pair  of  range  (R) 
u  -  variable  representing  the  Fourier  pair  of  velocity 
dfd  -  Doppler  frequency  resolution 
<5v  -  velocity  resolution 
SR  -  range  resolution 

Typically,  Doppler  processing  in  radar  consists  of  creating  range-Doppler  (R/D) 
maps — these  are  2-D  reflectivity  maps  of  the  scene  under  investigation  by  the  radar 
sensor,  plotted  as  a  function  of  the  range  and  radial  velocity  (or,  equivalently, 
Doppler  frequency)  dimensions.  The  creation  of  the  R/D  map  can  be  seen  as  the 
application  of  a  matched  filter  along  these  2  dimensions.  To  determine  the  equation 
governing  this  matched  filter,  we  need  to  derive  the  radar  response  of  a  point  target. 
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Assume  that  a  point  target  is  moving  with  constant  radial  velocity  v  toward  the 
radar  during  a  CPI.  Then,  its  range  varies  as  r(t)  =  R-  vt ,  where  R  is  the  range  in 

T  T 

the  middle  of  the  CPI  and  t  mns  from  — to  —  .  Also,  assume  that  the 

2  2 

transmitted  waveform  consists  of  a  train  of  N  short  UWB  impulses  with  complex 
envelope  p(t)  and  carrier  frequency /o,  with  a  repetition  interval  TV: 


N- 1 

-4t)  =  X  P(’  ~  nTr  )exp(./2/Z/V)  . 


71=0 


(1) 


The  signal  received  from  the  point  target  considered  previously  is  proportional  to  a 
delayed  version  of  the  transmitted  signal,  up  to  a  complex  constant  factor.  We  write 
the  received  signal  as 


r  PT 
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Note  that  throughout  this  note,  we  work  with  complex  analytic  signals,  meaning  we 
consider  both  the  in-phase  and  quadrature  (I-Q)  signal  components.  Upon  down- 
conversion  to  baseband,  we  eliminate  the  factor  exp(y'2;zf0t)  and  obtain 


N- 1 
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At  this  point,  we  need  to  mention  that  this  theoretical  derivation  is  valid  for  any 
type  of  UWB  waveforms,  not  only  short  impulses  as  considered  in  Eq.  1.  When 
pulse  compression  waveforms  are  used  for  this  purpose,  pit)  represents  the  pulse 

envelope  obtained  after  compressing  the  received  signal. 

The  signal  z(t )  is  sampled  in  the  time  variable,  and  its  samples  z(tk  )  are  rearranged 
in  a  2-D  array  indexed  by  the  discrete  variables  tfm  =  tk  -  nTr  (fast  time)  and 
tsn  -  nTr  (slow  time;  note  that  the  radar  waveform  is  already  discretized  along  this 
dimension  by  its  pulsed  nature).  When  denoting  the  new  2-D  signal  s(tjm , tm ) 

returned  by  a  point  target  at  range  R  and  moving  with  velocity  v,  we  add  the 
arguments  R  and  v  to  its  description: 
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SPT  fin  ’F„  ’R’V) 
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-R  +  —^Lvt„, 


with  m  =  0, ... M  - 1  and  n  =  — — , ...  —  - 1 .  In  the  last  equation,  we  neglected  tfm 

whenever  it  appeared  together  with  tsn  inside  the  same  factor,  since  tfin  « tsn , 
which  amounts  to  saying  that  the  pulse  width  of  p(t)  is  assumed  much  smaller  than 
the  repetition  interval  TV.  Note  the  strong  similarity  between  this  equation  and 
Eq.  11  in  our  previous  note,1  where  a  Gaussian  envelope  was  chosen  for  the 
transmitted  pulse.  Let  us  take  the  discrete  Fourier  transform  (DFT)  of  sPT  along 
the  tfm  dimension,  whose  Fourier  variable  pair  is  the  discrete  frequency  fm: 


SPT(fm,tsn,R,v)=P(fm)exp 


r  (■ 

i 

v  v 


M/„,  +  /(»  )R  |  4 +  /o  )vts„  ^ 


J) 


(5) 


=  P{fm)Q{fm,tSn’R’V) 


where  fm  =m  — ,  with  m  =  - 1,  and  P(f)  is  the  spectrum  of  pit) . 


P(f )  is  typically  a  real-valued  function;  this  is  the  case  when  p(t)  is  a  symmetric 
pulse  envelope  as  well  as  when  p[t)  is  obtained  by  compression  of  the  original 
transmitted  waveform. 


The  transfer  function  of  the  matched  filter  is  designed  to  compensate  the  phase  of 
the  response  of  a  point  target  at  range  R  and  moving  with  velocity  v.  Therefore, 
assuming  P(  f )  is  real,  the  matched  filter  transfer  function  is  Q  ( fm ,  tsn  .R.v).  The 

R/D  map,  which  represents  the  output  of  the  matched  filter  for  the  pair  of 
coordinates  (R,  v)  and  the  input  S(fm,tsn )  (the  response  of  an  arbitrary  scene),  can 

be  expressed  as 


M-lN-l 


RDM  (R,  v)  =  X  S(fm„ pT  v) 

+  /»  )R  4"(  /',  +  /(,  R, 


m= 0  n=0 
M -IN -l 


TJHS(fnrtsn)eXP 


m= 0 n= 0 


(Jl 


(6) 


J) 
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Before  we  take  a  more  in-depth  look  at  Eq.  6,  we  compare  this  formula  to  that 
obtained  by  conventional  Doppler  processing  (as  considered  in  our  previous 
work1).  Starting  again  with  the  S(fm,tsn )  2-D  signal,  the  R/D  map  in  conventional 

Doppler  processing  is  obtained  according  to 


M-lN-l 


RDMc(R,v)=Y,I,S(f„.l„)^P 


m= 0 «=0 


AnfmR  Wovtsl 


\  c 


(7) 


J) 


Equation  7  can  be  clearly  interpreted  as  a  2-D  DFT4  of  the  S(fm,tsn )  signal. 

(Although  R  and  v  are  formally  written  as  continuous  variables  in  our  equations, 
they  always  take  on  discrete  values  in  a  practical  implementation.)  Thus,  the  sum 
over  m  amounts  to  a  DFT  along  the  frequency-range  dimension,  while  the  sum  over 
n  amounts  to  a  DFT  along  the  slow  time-velocity  dimension.  The  additional  phase 


factor  exp 


jWqR 


appearing  in  Eq.  6  is  unimportant  and  was  simply  included  to 


make  the  formula  more  symmetric.  However,  the  factor  exp 


•  f mts 


J- 


,  which 


is  missing  from  Eq.  7,  is  critical  for  Doppler  processing  with  UWB  waveforms.  Its 
omission  from  Eq.  7  means  that  the  matched  filter  fails  to  compensate  for  the 
corresponding  factor  present  in  Eq.  5  and  explains  the  R/D  map  artifacts  for  the 
method  described  in  our  previous  work.1  Notice  that  neglecting  this  phase  factor  is 
well  justified  in  the  case  of  narrowband  signals,  where  fm  «  f0 ;  however,  this 

approximation  is  not  valid  for  UWB  waveforms,  where  fin  and/o  are  of  the  same 
order  of  magnitude. 

It  is  also  interesting  to  rewrite  the  formula  in  Eq.  6  as  a  function  of  the  delay  ( z) 

2  R 

and  Doppler  frequency  (fd)  variables,  which  are  connected  to  R  and  v  by  r  =  — 

c 

2  f 

and  f d  =  — —  v ,  respectively. 


M— IV— 1 


f  f 

RDM(t, f d )  =  XZ5(/'*,^)exp  ./'2^  (fm  +fo)T 

771=0  72=0  \ 


f  f  ^ 

1  +  IlL 


fo 


fdU 


J  ) 


■  (8) 


Going  back  to  Eq.  6,  we  introduce  the  new  variables  pm  =  +  —  1  and 


2(/m  +/oK 


and  rewrite  the  R/D  map  formula  as 
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(9) 


M-lN-1 

RDM  (R,v)=YjYjS  \Pm  ’  Um,n  )^v{i^{p,nR  +  «»,/))  • 

m=0  n=0 

The  last  equation  is  reminiscent  of  a  2-D  DFT  of  the  signal  s(pm,um  n),  involving 
the  Fourier  variable  pairs  (p,R)  and  (u,  v) .  However,  when  the  frequency  and  slow 
time  variables  are  sampled  on  a  regular  rectangular  grid  (at  equal  intervals  in  both 
dimensions),  the  samples  of  p  and  u  cover  a  wedge-shaped  domain  in  the  p  -  u 
plane,  as  illustrated  in  Fig.  1.  As  a  consequence,  we  cannot  apply  efficient 
algorithms,  such  as  the  2-D  fast  Fourier  transform  (FFT),  to  directly  compute  the 
expression  in  Eq.  9;  instead,  we  need  to  apply  a  brute-force  approach  to  evaluate 
Eq.  9.  A  further  discussion  of  the  computational  aspects  of  this  procedure  is  offered 
in  Section  4.1. 


Frequency  (GHz)  ,,  (i/m) 

(a)  (b) 


Fig.  1  Sample  points  used  in  creating  R/D  maps  by  the  matched  filter  method,  illustrating 
the  radar  data  support  in  (a)  the  frequency-slow  time  domain  and  (b)  the  p-u  domain.  The 
samples  correspond  to  the  numerical  example  in  Section  3,  down-sampled  for  figure  clarity. 


A  computationally  efficient  processing  scheme  for  the  implementation  of  Eq.  9, 
which  takes  advantage  of  the  2-D  FFT  algorithm,  can  be  obtained  using  the 
Keystone  transform.3  This  consists  of  remapping  the  radar  data  samples  from  their 
original  format  (rectangular  in  the/-  ts  plane)  onto  a  new  grid  such  that  the  resulting 
samples  in  the  p-u  plane  have  a  rectangular  support.  The  formula  describing  the 
coordinate  transformation  in  the/-  ts  plane  is 


t>L±A 

fo  2  /„ 


(10) 


When  replotted  in  the  /  -  ts '  coordinates,  the  radar  data  samples  are  arranged  in  the 
pattern  shown  in  Fig.  2a.  Correspondingly,  the  new  data  support  in  the  p  -  u  plane 


Approved  for  public  release;  distribution  unlimited. 


6 


takes  on  a  rectangular  shape  (Fig.  2b),  which  allows  the  application  of  the  2-D  FFT 
to  compute  the  expression  in  Eq.  9.  An  intermediate  step  consisting  of  interpolating 
the  radar  data  from  the/-  ts  to  the  /  -  ts '  coordinates  is  required  by  this  method. 


(a)  (b) 


Fig.  2  Sample  points  used  in  creating  R/D  maps  by  the  Keystone  transform  method, 
illustrating  the  radar  data  support  in  (a)  the  frequency-slow  time  domain  and  (b)  the  p-u 
domain.  The  samples  correspond  to  the  numerical  example  in  Section  3,  down-sampled  for 
figure  clarity. 

With  regard  to  the  resolution  of  this  radar  data-processing  approach,  the  same 
formulas  as  for  conventional,  narrowband  Doppler  processing  apply.  Thus,  the 
range,  Doppler  frequency,  and  velocity  resolutions  are,  respectively: 


SR  =  —  , 

2  B 

(11a) 

Sf\i  =  T  ’  and 

(lib) 

*  CPI 

<5v  -  C  • 

(He) 

m/o 


All  these  resolution  formulas  are  approximate  and  do  not  account  for  possible 
windowing  of  the  radar  data  in  the  frequency  and  slow  time  domains.  The  formulas 
in  Eqs.  1  lb  and  lie  assume  that  the  slow  time  interval  covered  by  the  radar  samples 
at/o  is  Tcpi.  Importantly,  the  resolution  in  both  dimensions  (R  and  v)  is  uniform 
over  the  entire  R/D  space,  unlike  the  processing  scheme  in  our  previous  work.1 

Before  concluding  this  section,  we  should  mention  that  the  Doppler  processing  with 
UWB  pulses  described  in  this  note  is  just  another  example  of  multidimensional 
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signal  processing  encountered  in  radar  applications.  Particularly,  we  can  draw  a 
parallel  with  the  formation  of  wideband,  widebeam  synthetic  aperture  radar  (SAR) 
images,  where  a  similar  coupling  between  the  2  dimensions  of  the  radar  data  space 
(frequency  and  aperture  position)  requires  the  application  of  a  2-D  matched  filter. 
The  conventional  Doppler  processing  is  analogous  to  the  Doppler  beam 
sharpening5  method  in  SAR  imaging,  which  is  known  to  produce  significant 
artifacts  when  used  in  conjunction  with  UWB  signals  and  wide  angle  apertures.  The 
brute-force  matched  filter  for  R/D  maps  described  by  Eq.  6  corresponds  to  the 
matched  filter  method  for  SAR  image  formation,  which  is  the  frequency-domain 
version  of  the  popular  backprojection  algorithm.5  On  the  other  hand,  the  Keystone 
transform  method  corresponds  to  the  range  migration  algorithm  for  SAR  image 
formation,6  which  uses  the  Stolt  interpolation  to  remap  the  data  samples  from  the 
original  wedge-shaped  support  in  the  kx  -  ky  plane  to  a  rectangular  grid,  followed 
by  the  application  of  a  2-D  FFT  (note  that  the  Stolt  transform  formula  differs  from 
the  Keystone  transform  due  to  different  sensing  geometries  and  radar  data 
dimensions,  but  the  principle  of  the  2  methods  is  very  similar). 


3.  Numerical  Example 


We  present  a  numerical  example  to  illustrate  the  processing  scheme  developed  in 
the  previous  section  by  a  simple  simulation.  Consider  2  point  targets  moving  with 
radial  velocities  vi  =  0.4  m/s  and  V2  =  -  0.6  m/s  (positive  v  indicates  a  target  moving 
toward  the  radar),  with  initial  ranges  R0l  =  4.05  m  and  RQ2  =3.3  m  .  The  radar 
transmits  a  train  of  modulated  Gaussian  pulses1  with  the  envelope 


p(t)  =  exp 


f  2\ 
t 


\  T0  ) 


.  The  pulse  width  is  r0  =  500  ps  (which  corresponds  to  a 


bandwidth  of  approximately  2  GHz)  and  the  carrier  frequency  is  /0  =  2.5  GHz .  An 


illustration  of  the  individual  pulse  in  the  time  and  frequency  domains  is  shown  in 
Fig.  3. 


The  pulse  repetition  interval  is  T,  =  25  ms,  which  corresponds  to  a  pulse  repetition 
frequency  (PRF)  of  40  Hz.  The  CPI  has  a  duration  Tcpi  =  1.2  s  and  involves  the 
coherent  processing  of  N  =  48  pulses.  The  total  observation  interval  is  3  s.  A 
representation  of  the  radar  data  returned  by  the  targets  as  a  function  of  range  and 
slow  time  is  shown  in  Fig.  4.  In  the  subsequent  processing  examples,  we  create  a 
single  R/D  map  corresponding  to  a  CPI  centered  at  t  =  1.5  s,  with  the  targets  at 
ranges  Rx  =  R()[  —vlt  =  3.45  m  and  R2  =  R02  —v2t  =  4.2m,  in  the  middle  of  the 

CPI.  The  unambiguous  velocity  interval  covered  by  this  processing  scheme  runs 
from  -1 .2  to  1.2  m/s. 
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Fig.  3  Representation  of  the  modulated  Gaussian  UWB  impulse  used  as  excitation  in  the 
numerical  example,  in  (a)  the  time  domain  and  (b)  the  frequency  domain.  The  red  dotted  lines 
represent  the  pulse  envelope. 


Range  (m) 


Fig.  4  Representation  of  the  radar-received  data  in  the  range-slow  time  domain,  showing 
the  2  targets  considered  in  the  numerical  example 

Each  received  waveform  within  a  repetition  interval  is  discretized  into  M  =  300  fast 
time  samples,  using  a  time  step  At  =  167  ps  .  The  total  fast  time  interval  covered  by 
these  samples  is  50  ns,  which  corresponds  to  a  target  range  running  from  0  to  7.5  m. 
The  frequency-domain  increment  is  20  MHz,  covering  6  GHz  of  bandwidth  (from 
1  to  7  GHz).  However,  since  the  processing  scheme  uses  I-Q  data  in  the  fast  time 
domain,  only  the  first  half  of  the  spectrum  (from  1  to  4  GHz)  contains  nonzero  data 
samples.  Note  that  the  frequency  offset  (fm  in  Eq.  6)  runs  from  -1.5  to  4.5  GHz, 
although  only  data  corresponding  to  frequency  offsets  from  -1.5  to  1.5  GHz 
contribute  to  the  sum  over  m. 
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To  mitigate  the  sidelobes  in  the  R/D  maps,  we  use  the  original  Gaussian  frequency- 
domain  window  of  the  transmitted  pulse  (for  range  sidelobe  suppression)  and  a 
Hanning  slow  time  domain  window  (for  Doppler  sidelobe  suppression). 


In  Figs.  5-7  we  present  the  R/D  map  obtained  for  the  CPI  centered  at  1.5  s  using  3 
different  processing  schemes:  1)  conventional  Doppler  processing,  corresponding 
to  Eq.  7  (Fig.  5);  2)  the  direct  implementation  of  the  matched  filter  described  by 
Eq.  6  (Fig.  6);  and  3)  the  matched  filter  implemented  via  the  Keystone  transform 
method  (Fig.  7).  The  target  images  in  Fig.  5  illustrate  the  main  issues  with 
conventional  Doppler  processing  applied  to  UWB  radar:  poor  resolution  in  both 
range  and  Doppler,  which  degrades  even  further  with  increasing  target  velocity.  On 
the  other  hand,  Figs.  6  and  7  demonstrate  that  careful  implementation  of  the  R/D 
matched  filter  for  UWB  radar  can  reach  the  same  Doppler  resolution  limit  as  a 
narrowband  radar  system,  independent  of  the  target  velocity.  The  only  difference 
between  the  target  images  in  Figs.  6  and  7  consists  of  the  sidelobe  directions,  which 
cannot  be  distinguished  in  the  R/D  maps  due  to  the  limited  dynamic  range. 


15 

10 
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-5 


-10 


-15 


-20 


Fig.  5  R/D  map  of  the  2  targets  considered  in  the  numerical  example,  obtained  by 
conventional  Doppler  processing,  for  a  1.2-s  CPI  centered  at  t  =  1.5  s 
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1 


20 


Fig.  6  R/D  map  of  the  2  targets  considered  in  the  numerical  example,  obtained  by  the 
matched  filter  method,  for  a  1.2-s  CPI  centered  at  t  =  1.5  s 


0  1  2  3  4  5  6  7 

Range  (m) 


Fig.  7  R/D  map  of  the  2  targets  considered  in  the  numerical  example,  obtained  by  the 
Keystone  transform  method,  for  a  1.2-s  CPI  centered  at  t  =  1.5  s.  Cubic  spline  interpolation 
was  used  in  the  implementation  of  the  Keystone  transform. 
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4.  Discussion 


4.1  Computational  Aspects 

In  this  section,  we  analyze  the  computational  aspects  of  the  various  UWB  Doppler 
processing  schemes.  As  suggested  by  a  comparison  of  Eqs.  6  and  7,  the 
implementation  of  the  exact  matched  filter  should  come  at  an  increased  complexity 
in  computational  terms,  as  compared  to  the  conventional  Doppler  processing 
formula. 

We  assume  that  the  radar  data  are  already  organized  in  a  2-D  array  of  size  M  by  N , 
S(fm,tsn ) >  with  the  2  dimensions  representing  frequency  and  slow  time.  Note  that 
some  UWB  waveform  implementations,  such  as  stepped  frequency  or  linear 
frequency  modulation  combined  with  stretch  processing,7  generate  radar  data 
directly  in  the  frequency  domain.  For  other  types  of  waveforms,  such  as  UWB 
impulses  or  pulse  compression  waveforms,  the  radar  data  are  generated  in  the  time 
domain  and  need  to  be  transformed  to  the  frequency  domain  via  FFTs. 

The  conventional  Doppler  processing  described  by  Eq.  7  is  efficiently  implemented 
as  a  2-D  FFT.  The  typical  form  of  this  algorithm  is  the  row-column  decomposition4 
and  requires  0(MN  log2  MN)  floating-point  operations  (the  symbol  O  stands  for 
“on  the  order  of ’)■  When  the  radar  data  are  directly  available  in  the  time  domain, 
there  is  no  need  to  convert  to  frequency  domain  in  this  Doppler  processing  scheme. 
Instead,  the  R/D  map  can  be  simply  obtained  by  1-D  FFTs  along  the  slow  time 
dimension  for  each  range  bin,  with  a  computational  load  of  0(MN  log2  N) . 
However,  one  should  keep  in  mind  that  the  difference  between  0{MN  log2  MN) 
and  0(MN  log 2  N)  is  smaller  than  it  appears  at  first  sight:  when  M  =  N,  the 
difference  is  a  factor  of  2. 

For  the  direct  (or  brute-force)  matched  filter  implementation  described  by  Eq.  6, 
we  cannot  apply  FFT-like  algorithms  due  to  the  shape  of  the  support  in  the  p  -  u 
domain.  Instead,  one  needs  to  compute  the  double  sum  in  Eq.  6,  pixel-by-pixel. 
Since  there  are  MN  pixels  in  an  R/D  map,  and  each  pixel  requires  0(MN) 

operations,  the  total  computational  load  is  o[m  2/V2),  which  is  significantly  larger 
than  that  required  by  the  2-D  FFT  algorithm.  Note  that  the  same  load  would  be 
obtained  if  we  implemented  the  2-D  DFT  in  Eq.  7  directly  by  double  summation, 
instead  of  taking  advantage  of  the  efficiencies  afforded  by  the  FFT  algorithm.  When 
the  original  radar  data  are  in  the  time  domain,  additional  0(MN  log2  M )  operations 

are  required  to  convert  them  to  the  frequency  domain  via  FFTs. 
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Another  option  for  the  direct  matched  filter  implementation  is  to  organize  the  data 
as  a  single-column  vector  of  length  MN,  using  the  Kronecker  product.7  In  that  case, 
the  R/D  map  is  also  organized  as  a  vector  of  length  MN,  while  the  matched  filter 
can  be  represented  as  an  MN  by  MN  matrix.  The  R/D  map  vector  is  then  obtained 
by  the  matrix-vector  multiplication  of  the  matched  filter  with  the  radar  data  vector. 
The  matrix- vector  multiplication  dictates  a  computational  load  of  o[m  2N2 ),  same 
as  the  method  described  in  the  previous  paragraph.  However,  this  method  is  not 
particularly  efficient  in  terms  of  computing  resource  usage  due  to  the  large  storage 
size  of  the  matched  filter  matrix,  as  well  as  the  time  required  to  compute  its 
elements  (as  an  example,  for  the  parameters  of  the  simulation  in  Section  3,  the 
matrix  would  have  a  size  of  1.66  GB).  This  approach  could  only  be  justified  in  a 
scenario  where  we  must  repeatedly  apply  the  same  matched  filter  to  a  stream  of 
radar  data,  such  as  in  creating  a  joint  range-time-frequency  representation 
(JRTFR)1  where  successive  R/D  maps  obtained  at  different  time  intervals  are 
stacked  into  a  3-D  data  structure. 

The  matched  filter  implemented  via  the  Keystone  transform  is  more  efficient  than 
the  brute-force  method,  since  it  uses  the  2-D  FFT  algorithm.  As  in  the  conventional 
Doppler  processing,  creating  the  R/D  map  requires  only  0(MN\og2MN ) 

operations,  which  represents  significant  savings  versus  the  o[m2N2)  operations 
required  by  the  brute-force  implementation.  However,  now  we  have  to  take  into 
account  the  extra  operations  involved  in  the  interpolation  required  by  remapping 
the  data  samples.  The  additional  computational  load  depends  on  the  interpolation 
method  employed;  this  can  vary  from  0(MN)  for  linear  interpolation  (which  is 
simple  but  generally  inadequate  for  this  purpose)  to  0(MN  log2  N)  for  Fourier-  (or 
sine-)  type  interpolation.  In  the  latter  case,  the  data  interpolation  step  approaches 
the  load  of  the  2-D  FFT  itself,  which  means  that  it  may  represent  a  large  part  of  the 
total  computational  time  required  by  this  method. 

4.2  Range  and  Doppler  Ambiguities 

So  far,  we  have  not  discussed  the  issue  of  ambiguities  created  by  the  range-Doppler 
processing  of  radar  data.  Since  this  is  an  important  problem  in  the  implementation 
of  radar  waveforms  and  processing  schemes,  this  section  analyzes  it  in  detail.  In 
particular,  as  suggested  by  Le  Chevalier,2  the  UWB  Doppler  processing  presented 
in  this  note  may  significantly  reduce  the  Doppler  ambiguities  as  compared  to  the 
conventional  processing  method.  The  derivation  of  the  AF  presented  in  this  section 
confirms  this  fact  in  a  quantitative  way. 

We  start  by  analyzing  the  conventional  (or  narrowband,  hence  the  subscript  N)  AF, 
as  defined  by  Rihaczek.8 
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(12) 


Xn  (a  fa )  =  {  (/  -  r)exp(-  j2nfdt)dt . 

-oo 

When  the  radar  waveform  x(t)  consists  of  a  train  of  baseband  pulses,  described  as 

N- 1 

x(t )  =  "y  p{t  -  nTr ) ,  this  function  becomes 

71=0 


N-l  » 

XN{TJd)=Yj \y)p\t-nTr  -T)exp(-  j2rfdt)dt .  (13) 

n=  0  — oo 

A  straightforward  interpretation  of  this  equation  is  that  the  AF  represents  the  R/D 
map  of  a  stationary  (v  =  0)  point  target  placed  at  range  R  =  0,  obtained  by  the 
conventional  Doppler  processing  method.  Equation  13  involves  pulse  compression, 
which  corresponds  to  correlating  the  received  signal  with  the  transmitted  pulse  p(t) . 

However,  as  discussed  in  Section  2,  our  method  of  creating  R/D  maps  omitted  the 
pulse  compression  step  by  considering  p(t)  as  either  a  UWB  impulse  or  the  output 

of  the  pulse  compression  operation  on  the  original  transmitted  waveform.  A  modified 
AF  that  more  closely  resembles  the  conventional  Doppler  processing  scheme 
described  in  Section  2  is  the  following  (note  the  superscript  S  for  “simplified”): 

N- 1  ® 

Xn  (a  fd  )  =  X  1  -  nTr  -  r)exp(-  jl7tfdt)dt 

n=0  _oo 

=  exp(-  jlrf dz)^x{j  +  nTr  )exp(-  j2nfdnTr ) .  (14) 

72=0 

The  only  change  in  Eq.  14  as  compared  to  Eq.  13  is  replacing  p(t)  inside  the 
integral  with  a  delta  function  S(t) .  The  second  expression  in  Eq.  14  can  be  further 
analyzed  by  the  methods  developed  by  Rihaczek.8  Letting  fd  =  0,  we  obtain  the 

N- 1 

expression  ^x(t  +  nTr),  which  describes  N  replicas  of  the  original  transmitted 

signal  along  Che  range  axis,  spaeed  apart  by  f ;  these  replieas  generate  «he  range 

ambiguities  in  the  R/D  map.  These  ambiguities  are  characteristic  to  any  radar 
waveforms  involving  a  periodic  train  of  pulses. 

In  this  work,  we  are  more  focused  on  analyzing  the  Doppler  ambiguities  obtained 
under  the  various  processing  schemes.  To  investigate  these,  we  let  r=  0  in  Eq.  14 

N- 1 

and  write  v(t)  explicitly  as  ^  p(t  - kTr ): 

*=0 
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(15) 


|  Zn  (0,  fd  )|  =  XZ  p((U  ~  k)T>  )eXP(“  nTr  )  • 


n=  0  k=  0 


Since  the  duration  of  p(t)  is  much  smaller  than  7V,  the  quantity  under  the  double 
sum  is  nonzero  only  when  n  =  k\ 


p((n-k)T) 


j p(o)  for  n  =  k 
[0  for  n  ^  k 


(16) 


Consequently,  we  obtain 


Xn  (O’  fd\  =  Yj  M°)exP(-  j24dnTr )  =  p{ 0) 


n=0 


sin  (rfdNTr) 
sin  (tfdTr) 


(17) 


The  expression  on  the  right-hand  side  of  Eq.  17  is  a  form  of  the  Dirichlet  (or  digital 
sine)  function,  frequently  encountered  in  digital  signal  processing  involving  DFTs. 
A  plot  of  ZN^’fd)  (cut  through  the  AF  along  the  Doppler  frequency  axis), 


obtained  for  the  simulation  parameters  in  Section  3,  is  shown  in  Fig.  8.  Notice  the 
equal  peaks  of  amplitude  Np(o) ,  spaced  apart  by  —  ;  these  peaks  (also  known  as 


grating  lobes)  generate  the  Doppler  ambiguities  in  the  R/D  maps  obtained  by 
conventional  Doppler  processing. 


Fig.  8  Cut  through  the  narrowband  AF  along  the  Doppler  frequency  axis,  demonstrating 
the  periodic  Doppler  ambiguities 
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While  the  AF  definition  in  Eq.  12  is  adequate  for  describing  Doppler  processing 
with  narrowband  waveforms,  it  does  not  capture  all  the  phenomenology  of  sensing 
a  moving  target  by  a  radar  system  operating  with  wideband  waveforms.  In  order  to 
account  for  the  latter,  we  introduce  the  so-called  “wideband  AF”,9  defined  as 


zw(T’fj)=  { 4ty 


-00  V 


f  f  ' 

t-T  +  —t 
fo  J 


exp(-  j2nfdt)dt . 


(18) 


N- 1 

For  a  train  of  baseband  pulses  x{t)  =  p(t  -  nTr ) ,  we  have 

n= 0 


N- 1  °° 


Xw 


4/4= X  J*4V 


f  r  \ 

t-nT„  -T  +  —t 

/c 


n—0  —co  \ 


exp(-  j2nf'dt)dt . 


(19) 


o  J 


Change  the  integration  variable  to  tf  =  t  -  nTr : 
Xw  (a  fd  ) 


V-1  ®  . 

=  X  1  XVf  +nTr)P 


(  f  r  \ 

i  +  Ll 


n= 0  . 
N-l 


fo 


fd 


A 


V  V  Jo  J 
iy  _x  00 

=  X  eXP(“  J2/fdnTr  )  J  4 /  +  )p 


-  r  +  —  nT 
fo  j 


n= 0 


tf  —  t  H — —nT 

’  fo  , 


ex  pi-jltfAf  +nTPlf 

exp(-  j2rfdt  f  ]dtf  ■  (2°) 


fd 


Similarly  to  the  narrowband  AF  analysis,  we  introduce  a  simplified  version  of  the 
wideband  AF,  which  closely  reflects  the  “wideband”  matched  filter  processing  used 
in  creating  R/D  maps  with  UWB  waveforms  (Eq.  6).  This  new  AF  is  obtained  by 
replacing  p(t)  inside  the  integral  with  a  delta  function  S(t) : 


Xw  (A  fd  ) 


N-l 


n= 0 


=  X  exp(“  j27f d  n  Tr  )  |  4. /  +nTr)8  tf~T  +  'ynTr  ]  eXP  (“  4 

-00  V  J  0  / 


N-l 


exp(-  x\  r  +  nTr  -  J-  nT  exp(-  j2tfdnTr ) 

Jo  ) 


n= 0 


.  (21) 


N-l 


After  replacing  x{t)  by  ^  p(t  -  kTr ) ,  we  obtain 


k= 0 


N-l  N-l 


n= 0  k= 0 


| Xw  (a  fd  )|  =  X  X  p\  T  +  (n  -  k  )Tr  -JfnT,  exp(-  j2rfdnTr )  •  (22) 

Jo  J 
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By  a  similar  argument  as  the  one  used  in  Eq.  16  and  noting  that  fd  «  f0 ,  we  can 
write 


r  +  (n  -  k  )Tr 


f 

-—nT 

fo  J 


f p(r)  for  n  =  k 
[0  for  77  ^  k 


(23) 


Consequently, 


TV  1 


xw{rJd)  =  YjP 


77=0 


r  -—nT„ 


/< 


exp(-  j2nfdnTr). 


(24) 


o  J 


Now  let  us  write  p 


fd 


as  an  inverse  Fourier  transform,  using  rand  fm  as 


r  -  —  nT 
v  fo  y 

the  Fourier  variable  pair  (note  that/,,,  are  discrete  frequency  samples): 


f 


fd 


\ 


t  -  —  nT 
fo  j 


M- 1 


r 


m= 0 


f 


=  Z  P(f»<  )eXP|  J27r\  (fm  +  fo  y  -  fm  ^fnT, 

Jo  J) 


W 


v  v 


(25) 


Then,  the  wideband  AF  can  be  written  as 


M- 1  N-\ 


xsw(*Jd)  =  XXp(/,Jexp 


m= 0  72=0 


r  r 

jin 

\  v 


(/m+ZoV- 


1  +  4^ 

v  /o; 


A3 


(26) 


We  can  interpret  this  form  of  the  wideband  AF  as  the  output  of  the  “wideband” 
matched  filter  (as  in  Eq.  6)  for  a  stationary  point  target  placed  at  R  =  0.  Thus,  the 
response  of  this  target  can  be  written  as 


Sfr(/..<».0.0)=  />(/„). 


(27) 


The  R/D  map  generated  by  this  response,  as  a  function  of  rand  fd  ,  is  (see  Eq.  8) 


M-lN-l 


f  / 

RDM (t,  fd  )  =  Z  Z  P(Z»  )exP  3'2yr  (/„,  +  /0  )r 

V  v 


772=0  72=0 


1  +  4^ 

V  fo  J 


f dfr 


w 


J) 


(28) 


which  is  identical  to  the  expression  in  Eq.  26  given  the  fact  that  tsn  -  nTr . 

To  investigate  the  Doppler  ambiguities,  we  go  back  to  Eq.  24  and  take  a  cut  through 
the  AF  along  the  Doppler  frequency  axis  (or  r=  0): 
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cxp(-  j2nfdnTr). 


(29) 


N- 1 


xw(°’fd}  =  TjP 


n= 0 


f  ) 

—  nT 

fo  J 


It  is  clear  that  this  function  does  not  resemble  the  digital  sine  obtained  in  the 
narrowband  case.  In  fact,  when  we  plot  it  as  a  function  of fd  (again  for  the  simulation 
parameters  in  Section  3),  we  obtain  the  graph  in  Fig.  9.  Notice  that  the  peaks 

(grating  lobes)  present  at  fd  =  ± —  are  strongly  attenuated  with  respect  to  the  peak 
at  fd  =  0.  We  quantify  this  attenuation  by  computing  the  ratio 


Zw 

f  n 
l  Tr 

N-l 

I > 

n= 0 

f  \ 

n 

K  y 

|rf(0,o|  Np( 0) 

(30) 


Fig.  9  Cut  through  the  wideband  AF  along  the  Doppler  frequency  axis,  demonstrating  the 
attenuation  of  the  Doppler  ambiguities 


Since  p(t)  is  monotonically  decreasing  away  from  t  =  0,  it  follows  that  the  ratio  in 
Eq.  30  is  always  smaller  than  1 .  Moreover,  since  p(t)  has  a  narrow  support  (typical 


for  an  UWB  pulse),  we  expect  p 


\ 

n 

fo  y 


to  be  zero  for  n>  N{,  where  Ni  is  an 


integer  smaller  than  N.  Then  the  first  grating  lobe  attenuation,  or  the  ratio  in  Eq.  30, 

can  be  approximated  by  — - .  We  can  also  establish  an  approximate  relationship 

N 

between  Ni  and  the  pulse  width  of  p(t) ,  denoted  by  r0 :  —  =  r0  =  — ,  where  B  is 


fo 


B 


the  pulse  bandwidth. 
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The  condition  required  to  obtain  a  large  attenuation  of  the  first  grating  lobe  is 

^L«1 ;  by  replacing  N\  with  —  (or  the  inverse  of  the  fractional  bandwidth),  we 
N  B 


finally  obtain  N  » 


fo 

B 


This  is  identical  to  the  formula  indicated  by  Le  Chevalier.2 


Translated  in  words,  this  formula  says  that  the  number  of  pulses  in  a  CPI  must  be 
much  larger  than  the  inverse  of  the  fractional  bandwidth  to  achieve  strong 
attenuation  of  the  Doppler  grating  lobes  in  the  wideband  AF.  This  condition  is 
relatively  easy  to  satisfy  for  UWB  pulses,  which  by  definition  have  a  large 

f 

fractional  bandwidth  (or  small  —  ratio). 


To  clearly  demonstrate  the  effect  of  the  conventional  and  wideband  Doppler 
processing  schemes  on  the  Doppler  ambiguities,  we  revisit  the  numerical  example 
in  Section  3  by  extending  the  velocity  axis  in  the  R/D  maps  to  run  from  -2.4  to 
2.4  m/s.  Regardless  of  the  processing  scheme,  both  targets  create  an  alias  visible  in 
these  maps:  the  first  target  at  -2  m/s  and  the  second  target  at  1.8  m/s.  However,  the 
amplitude  of  these  ambiguities  differs  dramatically  between  the  conventional 
Doppler  processing  (Fig.  10),  where  they  are  as  large  as  the  true  targets,  and  the 
wideband  Doppler  processing  (Fig.  11),  where  they  are  strongly  attenuated. 


0 

□ 

0 

o 


0  1  2  3  4  5  6  7 

Range  (m) 


Fig.  10  R/D  map  of  the  2  targets  considered  in  the  numerical  example,  obtained  by 
conventional  Doppler  processing,  showing  an  extended  velocity  scale  to  demonstrate  the 
Doppler  ambiguities 
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0  1  2  3  4  5  6  7 

Range  (m) 


Fig.  11  R/D  map  of  the  2  targets  considered  in  the  numerical  example,  obtained  by  the 
matched  filter  method,  showing  an  extended  velocity  scale  to  demonstrate  the  Doppler 
ambiguities 

The  importance  of  this  result  in  improving  radar  performance  can  be  appreciated  in 
the  context  of  2  traditional  tradeoffs  typically  encountered  in  radar  systems:  1)  the 
tradeoff  between  range  and  Doppler  resolutions8  and  2)  the  tradeoff  between  range 
and  Doppler  ambiguities.5  Thus,  the  shortcomings  of  the  conventional  Doppler 
processing  proposed  in  our  previous  work1  seemed  to  indicate  that  we  cannot 
simultaneously  increase  the  range  and  Doppler  resolutions  past  certain  limits  by 
choosing  both  large  signal  bandwidth  and  long  CPI  duration  due  to  the  issue  of 
target  range  migration  within  a  CPI.  However,  these  conflicting  requirements  are 
resolved  through  the  wideband  matched  filter  processing  scheme.  This  allows  us  to 
choose  a  long  CPI  for  good  Doppler  resolution,  without  limiting  the  signal 
bandwidth  (which  dictates  the  range  resolution). 

Regarding  the  tradeoff  between  range  and  Doppler  ambiguities,  this  is  dictated  by 
the  choice  of  PRF.  In  conventional  Doppler  processing,  avoiding  the  grating  lobes, 
which  occur  at  intervals  equal  to  PRF,  typically  requires  a  large  PRF.  On  the  other 

c 

hand,  avoiding  the  range  ambiguities,  which  occur  at  intervals  equal  to  9p^p » 

requires  a  small  PRF.  Several  complex  waveforms  and  signal  processing 
techniques  have  been  developed  to  mitigate  this  tradeoff5,8  (e.g.,  transmission  of 
trains  of  pulses  with  variable  PRF),  most  of  them  involving  additional  penalties  in 
system  performance.  Nevertheless,  the  processing  scheme  proposed  in  this  note 
resolves  this  tradeoff  by  naturally  suppressing  the  Doppler  grating  lobes  as 
demonstrated  in  this  section.  Note  that  this  effect  holds  only  for  the  case  of  UWB 
waveforms  (i.e.,  waveforms  with  large  fractional  bandwidth). 
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5.  Conclusions 


In  this  note,  we  presented  a  UWB  radar  Doppler  processing  scheme  that  overcomes 
the  limitation  reported  in  some  of  our  previous  work.  The  method  is  based  on  a 
careful  derivation  of  the  matched  filter  for  this  type  of  processing,  taking  into 
account  the  coupling  between  the  frequency  and  slow  time  dimensions  in  the  radar 
response  of  a  point  target.  This  coupling  leads  to  a  2-D  data  transformation,  similar 
to  that  employed  in  the  formation  of  wideband,  widebeam  SAR  images. 

Theoretical  and  numerical  comparisons  of  the  new  method  with  conventional 
Doppler  processing  clearly  demonstrate  the  superior  performance  of  the  wideband 
matched  filter  in  terms  of  both  range  and  Doppler  resolution,  and  at  the  same  time 
explains  why  the  old  approach  fails  when  UWB  waveforms  are  employed.  We  also 
present  an  efficient  implementation  of  the  new  technique,  based  on  the  Keystone 
transform,  followed  by  a  discussion  of  the  computational  complexity  of  all  these 
processing  schemes. 

The  detailed  derivation  of  the  wideband  AF  in  Section  4.2  allowed  us  to  confirm 
the  result  reported  by  Le  Chevalier2  regarding  the  mitigation  of  Doppler 
ambiguities  by  UWB  waveforms.  This  is  an  important  result  showing  that  it  is 
possible  to  resolve  some  traditional  tradeoffs  in  radar  system  performance  by  using 
this  type  of  waveform. 

Regarding  possible  applications  of  this  new  Doppler  processing  method,  one 
should  mention  the  combination  of  micro-Doppler  and  high-range  resolution 
analysis  by  means  of  the  JRTFR.  Given  the  good  resolution  in  all  dimensions 
(range,  Doppler,  and  time),  this  type  of  analysis  should  be  very  beneficial  to  target 
classification  problems.  An  excellent  example  of  this  processing  applied  to  the 
analysis  of  a  walking  human  target  is  presented  by  He  et  al.3  We  are  confident  that 
future  investigations  will  reveal  new  applications  of  the  wideband  Doppler 
processing  technique. 
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List  of  Symbols,  Abbreviations,  and  Acronyms 


1-D 

1 -dimensional 

2-D 

2-dimensional 

3-D 

3-dimensional 

AF 

ambiguity  function 

CPI 

coherent  processing  interval 

DFT 

discrete  Fourier  transform 

FFT 

fast  Fourier  transform 

I-Q 

in-phase  and  quadrature 

JRTFR 

joint  range-time-frequency  representation 

PRF 

pulse  repetition  frequency 

R/D 

range-Doppler 

SAR 

synthetic  aperture  radar 

UWB 

ultra-wideband 
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